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Revisiting the He n to H i ratio in the Intergalactic Medium 
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ABSTRACT 

We estimate the He n to H i column density ratio, rj = N(Heii)/N(Hi), in the intergalactic 
medium towards the high redshift (z e m = 2.885) bright quasar QSO HE 2347-4342 using 
Voigt-profile fitting of the H i transitions in the Lyman series and the He n Lyman-a' transition 
as observed by the FUSE satellite. In agreement with previous studies, we find that r\ > 50 
in most of the Lyman-ff forest except in four regions where it is much smaller (77 ~ 10 - 20) 
and therefore inconsistent with photo-ionization by the UV background flux. We detect O vi 
and C iv absorption lines associated with two of these regions (z a b s = 2.6346 and 2.6498). We 
show that if we constrain the fit of the H 1 and/or He 11 absorption profiles with the presence 
of metal components, we can accommodate 77 values in the range 15-100 in these systems 
assuming broadening is intermediate between pure thermal and pure turbulent. While simple 
photo-ionization models reproduce the observed A^(0 vi)/7V(C iv) ratio, they fail to produce 
low 77 values contrary to models with high temperature (i.e T > 10 5 K). The Doppler pa- 
rameters measured for different species suggest a multiphase nature of the absorbing regions. 
Therefore, if low j] values were to be confirmed, we would favor a multi -phase model in which 
most of the gas is at high temperature (> 10 5 K) but the metals and in particular C iv are due 
to lower temperature (~ few 10 4 K) photo-ionized gas. 
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The presence of metals in the H 1 Lyman-ff forest at optical depths 
TLya ^ L detected through C iv and O vi ab sorption lines seen 
in QSO spectra, is now well established (see [Songaila & Cowid 
1 1994 iBergeron et alJ[20o3 : ISimcoe et alj|2004h . Observations are 
consistent with an average carbon metallicity relative to solar 
of [C/H] ~ -2.8 with no sign of redshift evolution over the 
range 1.8 < z < 4.1 but a significant trend with over-densities 
dSchave et al.ll2003l ; lAguirre et al]|2008h . Given the expected low 
metallicities and the high ionization state of the gas, direct de- 
tection of metal absorption lines from underdense regions of the 
intergalactic medium (IGM) is beyond the scope of present day 
large te lescopes. Statistica l methods like pixe l analysis are used 
instead tellison et alJ|2qO0l:ISchave et alj|2003l ; I Aracil et alj|2004 
IScannapieco et alj|2006l : IPieri et al.l l2010h and show that metals 
must be present in the low-density regions. Even in regions where 
C iv absorption is detected directly, it is not clear however what is 
the main physical process that is maintaining the ionization state of 
the gas. In general, it is believed that photo-ionization keeps the gas 
ionized. However, it is probable that mechanical inputs from galac- 
tic winds can influence the ionization state of part of the IGM gas 
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through collisional ionization at least in the proximity of galaxies. 
Therefore, it is important to simultaneously study different species 
covering a wide range of ionization states to get a better understand- 
ing of the metal enrichment and the different ionizing mechanisms 
at play. 

Recent hydrodynamical simulations jDave et alj l200ll : 
iFang & Brva5l200ll : [Rang et alj|2005l iBouche et alj|2006l. (20071) 
suggest that the missing baryons at low redshift, z ~ - 0.5, and 
the missing metals at high redshift, z ~ 2.5, could reside in the 
warm-hot phase of the intergalactic medium (called WHIM) with 
T ~ 10 5 - 10 7 K). Highly ionized species of oxygen such as O vi, 
O vii and O vm can be useful probes of the WHIM. While the 
strongest transitions of the latter two species have rest-wavelengths 
in the soft X-ray range, the spectral doublet O vi/i/il032,1037 is 
seen in the near UV range and is therefore a useful probe of the 
gas at a temperature of ~3xl0 5 K, temperature at which O s+ /0 is 
maximum. 

It has been suggested that a large fraction of the conspicuous 
O vi phase seen to be associated to high-z damp ed Lyman-tt sys - 
tems may originate from collisionally ionized gas dFox et alj2007t) . 
However, photo-ionization can also maintain oxygen in a high ion- 
ization state and at relatively low temperature (T ~ a few 10 4 K, see 
Oppenheimer & Dave 2009). Actually a large fraction of the O vi 
absorption seen at z > 2.5 in quasars show Doppler broadening 
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consiste nt with photo-ionization in the vicinity of the QSOs (for 
example ISrianand & PetitiearJ boOO) but also in the diffuse inter- 
galactic medium (for example. iBergeron et alj|2002l : Isimcoe et al.l 
120041) . While it is expected that the intrinsic ionizing spectrum of 
QSOs is hard enough to maintain a high degree of ionization of 
oxygen in their vicinity, in the IGM, the hardness of the ionizing 
spectrum will depend on the intrinsic spectral shape of the ioniz- 
ing sources and the IGM opacity at the H 1 and He 11 Lyman Limit 
dHaardt & Madaull 19961: iFardal et al.lll998l) . 

An additional piece of information comes from QSO lines of 
sight transparent in the Lyman continuum [i.e a high-z QSO line 
of sight without any intervening Lyman limit system (LLS) block- 
ing the UV end of the spectrum]. It is then possible to observe 
the rest wavelength ranges of the H 1 and He 11 Lyman-o- forests 
and to compute the ratio of the Hen to the Hi optical depth (i.e 
the ri parameter). Th e bright QSO HE 2347-4342 at z em = 2.885 
( Reim ers et aljfl997h is one such targets that attracted a lot of at- 
tention in the past years. It was shown that the Hen opac ity is 
"patchy" in nature dReimers et~al]|l997l: [Smette et ai]|2002h and 
that t] decreases gradually from higher to lower reds hift possibly 
due to a change in the sl ope of the ionizing spectrum 1 Zhe ng et al.l 
120041) . Shull et al. (2004) discussed the small scale variations (over 
Az ss 10~ 3 ) of 77 and found an apparent correlation between high rj 
(less ionized He 11) and low H 1 column density. They ascribed these 
small scale rj variations to "local ionization effects" in the proxim- 
ity of QSOs located close to the li ne of sight and havin g spectral 
indices ranging from a s = to 3. IWorseck et al.l d2007l) reported 
the detection of 14 foreground QSOs in the field located close to 
the line of sight and could not find any convincing evidence for any 
transverse proximity effect from a decrease in the H 1 absorption, 
although they did claim that the local UV spectrum inferred in the 
vicinity of three foreground QSOs appeared harder than expected, 
which is an indication of a transverse proximity effect. In turn these 
fluctuations could be due to an appreciable contribution of thermal 
broadening to the velocity width of absorption lines at high N(H 1) 
dFechner & Reimersll2007l) . 

In this paper and after a description of the observations (Sec- 
tion 2), we use a different approach involving Voigt profile fitting 
analysis of the H 1 and He 11 absorption lines to measure 77 (Sec- 
tion 3). We then report new detections of O vi absorption associated 
with regions with low rj values (Section 4) and construct models of 
these regions (Section 5) before concluding in Section 6. 



2 OBSERVATIONS 

The optical spectrum of HE 2347-4342 (z cra = 2.885) used in 
this study was obtained wit h the VLT UV Echelle Spectrograph 
(UVES) dDekkeretal1l2000h mounted on the ESO Kueyen 8.2- 
m telescope at the Paranal observatory in the course of the 
ESO-VLT large prog ramme 'The Cosmic Evolution of the IGM' 
( Berg eron et atll2004h . HE 2347-4342 was observed through a 1- 
arcsec slit (with a typical seeing of 0.8 arcsec) for 12 h with cen- 
tral wavelengths adjusted to 3460 and 5800 A in the blue and red 
arms, respectively, using dichroic #1 and for another 14 h with cen- 
tral wavelengths at 4370 and 8600A in the blue and red arms with 
dichroic #2. The raw data were reduced u sing the latest version 
of the UVES pipeline (Ball ester et al.l l 2000) which is available as 
a dedicated context of the MIDAS data reduction software. The 
main function of the pipeline is to perform a precise inter-order 
background subtraction for science frames and master flat fields, 
and to apply an optimal extraction to retrieve the object signal, re- 



jecting cosmic ray impacts and performing sky subtraction at the 
same time. The reduction is checked step-by-step. Wavelengths 
are corrected to vacuum-heliocentric values and individual one- 
dimensional spectra are combined. Air-vacuum conversions and 
helioc entric corrections were don e using standard conversion equa- 
tions dEdlenl 1966lJstum P n1 1 980). Addition of individual exposures 
is performed by adjusting the flux in individual exposures to the 
same level and inverse variance weighting the signal in each pixel. 
Great care was taken in computing the error spectrum while com- 
bining the individual exposures. Our final error in each pixel is the 
quadratic sum of the weighted mean of errors in the different spec- 
tra and the scatter in the individual flux measurements. Errors in 
individual pixels obtained by this method are consistent with the 
rms dispersion around the best fitted continuum in regions free of 
absorption lines. The final combined spectrum covers the wave- 
length range of 3000 to 10, 000 A. A typical SNR-60 per pixel (of 
0.035 A) is achieved over the whole wavelength range of interest 
for a spectral resolution of R ~ 45, 000. The detailed quantitative 
des cription of data cal ibration is presented in lAracil et al. I d2004l) 
and lChand etafl d2004l) . 

We use the continuum normalized FUSE data provided by Dr. 
Zheng. The details o f the data reduction and continuum normaliza- 
tion can be found in|zheng et alj d2004h W The original data have 
typical resolution of R = 20, 000 and signal-to-noise ratio ~5 in 
the lo ng wavelength range (A > 1050 A). Following Zhe ng et"al] 
(2004), we have re-binned this data to 0.1 A, which leads to an 
effective resolution of R ~ 4000. We restrict ourselves to the wave- 
length range with SNR > 10. This corresponds to a redshift range 
2.58< z<2.70 or a velocity range of ~10, 000 km s around a cen- 
tral redshift of z = 2.6346 (see Fig.[T](. 



3 /V(He n)/2V(H 1) RATIO 

In this Section we concentrate on the column density ratio 
r] = /V(He n)/7V(H 1) over the redshift range 2.58< z < 2.70 where 
the FUSE data show relatively good signal to noise ratio. This range 
roughly corresponds to a relative velocity range of -4000 km s _1 to 
+ 5000 km s -1 around the strong O vi absorber seen at z a bs = 2.6346 
(see Fig.[T](. 

As a first step we fitted simultaneously the Lyman-ar to 
Lyman-y profiles when possible, e.g. when the Lyman-/? and/or 
Lyman-y lines are not blended with another intervening Lyman- 
a line. Then we compare the He 11 absorption profile with a model 
with the same components as the H 1 model, scaling the fitted H 1 
column densities by the parameter rj. We consider two alternatives: 
In the first case we use the same Doppler parameter for H 1 and 
He 11 (assuming turbulent broadening); in the second case we give 
the He 11 ^-parameter the value expected from thermal broadening 

1 We have obtained individual spectra reduced using Calfuse 3.2. 1 version 
from http://fuse.iap.fr/interface.php We combine LIF spectra after correct- 
ing for the background by demanding zero flux in the core of strong satu- 
rated absorptions in the wavelength range 1130-1185 A. When we follow 
the same continuum fitting and re-binning procedures, we find the new data 
follow the s tructures (both in wavelength and flux) as seen in the data of 
Zhen g et alj j2004h very well and fitting results are not changed. So, results 
presented in this paper will not change when one uses the new pipeline for 
the data reduction. Whereas this w ork was already co mpleted, new COS 
data on this object were reported by Shull et al. 12010). As the COS spec- 
trum is found to be consistent (see their Fig. 3) with the FUSE spectra used 
here, this has no consequence on the results of this paper. 
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Figure 1. H i and He n Lyman-a and H i Lyman-/? absorption profiles on a velocity scale with origin z = 2.6346. Different velocity ranges where i] is measured 
through x 1 minimization are indicated with horizontal lines and labelled by letters in alphabetical order from left to right. Vertical tick marks located above the 
absorption profiles show the positions of the individual Voigt profile components used to fit the Lyman-tr and Lyman-/? H i lines together. The best fit models 
for H i Lyman-/? and Lyman-a are overplotted to the data. The FUSE spectrum is shown in the top panel together with the best fitted profile obtained from 
scaling the H i column densities by the fitted rj parameter and assuming thermal (solid red line) or turbulent (dashed blue line) broadening for the He n lines. 



(i.e b(He n) = 0.5xfo(H i)). The best fitted values of 77 is obtained by 
X 2 minimization. While fitting the He 11 profiles we use a Gaussian 
convolving function to correctly represent the FUSE spectral res- 
olution. For Voigt -profile decomposit ion we have used the fitting 
code developed by Kh are et al. (1997). 

As the FUSE data are of much lower resolution and SNR than 
the UVES data, we cannot estimate 77 for individual H 1 compo- 
nents. Instead, we have singled out 15 small regions named as A, 
B, C etc., in Fig.[TJand we derive the best 77 value over each region. 
We wish to point out that the approach w e have taken here i s very 
different from previous studies. Indeed, IShull et alj (2004) used 
apparent optical d epth in Lyman-g only ( AOD) method, whereas 
iKriss et all d200ll) andlzheng et alj d2004h used Gaussian decom- 



position and Fechner & ReimersT i 20071) scale the whole H 1 spec- 
trum by r\ = 4xt Hc ii/thi to fit the He 11 data. In all these studies only 
H 1 Lyman-cf is used0. This is the use of the H 1 optical depths in all 
available Lyman series lines that allows us to discriminate between 
thermal and turbulent broadening. 

The best fitted Voigt profiles to the Lyman-ce and Lyman-/? 



2 limi ted amount of analysis of Lyman-/? have been done by Zh eng et"al] 
<2004h . 



absorption lines are shown in the bottom and middle panels of 
Fig- [3 The top panel shows the best fitted He 11 Lyman-o- line with 
the two assumptions on the Doppler parameter discussed above. 
The x 1 curves as a function of 77 for the different regions singled 
out in Fig. [TJ are shown in Fig. [2] The solid and dashed lines in 
these plots represent the cases of thermal and turbulent broadening 
respectively. In most cases the x 1 curve shows a clear minimum 
thereby allowing us to discriminate between the turbulent and ther- 
mal cases, and to derive the best fitted value of 77. Errors are esti- 
mated from the range of 77 values corresponding to Ax 2 = ± 1 around 
the minimum. There are regions, especially when the He 11 Lyman- 
a line is saturated, for which the x 2 curve flattens (e.g. regions E 
and G), we have only one-sided limit. In these cases we define the 
2cr lower limit of 77 as the value corresponding to ax 2 equal to^ 2 of 
the flat part of the curve plus four. The shapes of the x 2 curves are 
not symmetric which is a natural consequence of line saturation. 

It is clear from the Figure that, apart from region I, the x 2 
values are smaller in case of turbulent broadening and that minima 
are reached only in that case. In the case of thermal broadening, 
the x 2 curves seem to saturate to some asymptotic value probably 
because the observed He 11 profiles are too broad to be reproduced 
by the model. Thus the exercise presented here shows that the width 




Figure 2. x 1 resulting from the comparison of the He n absorption with a model profile obtained by scaling the H i column density with the parameter rj as a 
function of i] for the different regions defined in Fig.[T] The dashed blue and solid red curves are for the two extreme cases, respectively: b(He 11) = b(H i) (i.e 
turbulent broadening) or fo(He n) = 0.5xfo(H i) (i.e pure thermal broadening). 
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Figure 3. The r\ values for the 15 regions singled out in Fig. [T] are plotted 
against the neutral hydrogen column density. 

of He ii Lyman-o' lines are consistent with the fo-parameter derived 
from H i lines. 

If the gas is optically thin and photo-ionized by a UV back- 
ground dominated by QSOs, we would expect rj to be in the range 
40-400 depending on the exact spectral index and the IGM opac- 
ity. In t he case of self- shie lded optically thick gas, rj could be even 
higher dFardal et al.lll998h . Four regions (D, F, H and N) in Fig.|2] 
have j] < 40. These regions are associated with large H i column 
densities as can be seen on Fig. [3] where log rj is plotted against 
log Af(H i) as measured in the different regions. This correlation 
was already noted in earlier works. Fechner & Reimers (2007) ar- 
gued that this can be explained if the thermal broadening of lines 
are also important. 

In the following, we will use additional information on metal 
lines observed in the UVES spectrum to discuss further the ioniza- 
tion state of the gas in these regions. 



4 REGIONS WITH LOW 77 VALUES 

In the previous Section, we have shown that the N(He ii)/A^(H 1) 
ratio can be explained over most of the observed spectrum by ion- 
ization of the gas by the UV background except in four regions: 
D (z ab5 = 2.6346), F (z abs = 2.6498), H (z abs = 2.6624) and N (z abs = 
2.6910). The presence of O vi and C iv absorption in systems show- 
ing low values of rj may yield interesting clues about (i) the nature 
of the ionizing radiation, (ii) the effect of thermal/turbulent broad- 
ening and (iii) the possible mechanical feedback from winds. 

Regions D and F are associated with C iv and strong O vi ab- 
sorption lines. These are the only two Civ systems in the redshift 
range 2.58 < z < 2.70 (see top panels in Fig. [4j and we discuss 
them in detail below. 

For region H, the wavelength range where possible O wiA 103 1 
absorption is redshifted is strongly blended and only a possible 
weak line is present at the expected location of Ovi/11037 (see 



bottom-left hand-side panel in Fig. [4}. As no other metal line is 
detected at this redshift we are unable to confirm if this feature is 
indeed due to O vi absorption. 

For region N, while both O vi regions are strongly blended, the 
optical depth constraints are satisfied at two velocity positions (see 
Fig. U). However, the possible O vi/11031 feature is also consis- 
tent with being C m absorption at z abs ~2.8972. Similarly, the pos- 
sible O vi/11037 line is blended with Lyman-/? at z abs = 2.7306 and 
O vi/11031 at z abs = 2.7121. Hence, we cannot confirm the presence 
of Ovi absorption in this region. Note that in region N (i.e. z abs = 
2.6910) 77 is probably affected by transverse proximity effect from 
QSO J23495-4338 located at redshift z = 2.690±0.006, 15 arcmin 
away from the line of sight of interest JWorseck et alj2007h . 

Note that we detect Ovi absorption at z abs = 2.7121, 2.7356 and 
2.7456 as well. The He 11 opacity is highatz abs = 2.7121 and 2.7456 
which makes rj difficult to estimate. If we scale the Voigt-profile 
fits to the Hi absorption to reproduce the Hen profile, we find rj 
to be in the range 10-100 and >100 for, respectively, the systems 
at z abs = 2.7121 and 2.7456. The wavelength range in which the 
Hen absorption at z abs = 2.7356 is expected to be redshifted has 
been removed by Zhe ng et all J2004h because of the strong airglow 
lines so that we cannot estimate rj for this system. The system at 
z abs = 2.7356 is a known Lyman-limit system. A Voigt-profile fit to 
the Hi absorption gives log A'(Hi) = 16.50±0.28. 

4.1 System at z abs = 2.6498 

A velocity plot of high ionization metal lines and H 1 lines from this 
system is shown in the top-right panel of Fig. [4] Clearly the metal 
lines are off-centered with respect to the H 1 absorption. In addition, 
there is a velocity off-set of 2 to 10 km s -1 between the centroids of 
the O vi and C iv absorption profiles. Interestingly, all the shifts are 
in the same direction as would be expected in a flow ionized from 
the same side. The best fit of the profiles is obtained when we allow 
for C iv component redshifts to be independent of that of the O vi 
components (see Table[T](. 

Doppler parameters are larger for Ovi compared to Civ 
which supports neither pure thermal nor pure turbulent broaden- 
ing. The upper limits on the kinetic temperature of the gas mea- 
sured from the ^-parameters of O vi components is 1.4x10 s , 8xl0 4 , 
1.2x10 s and 4xl0 4 K respectively for components at -76.2, -48.6, 
-20.2 and +4.2 km s" 1 . Therefore within the allowed error in b- 
parameters, the O vi profile allows for the existence of high tem- 
perature (r > 10 5 K) at least in part of the associated gas. 

The top panel of Fig. [5] shows the apparent column densities 
of O vi (in blue) and C iv (in red) per unit velocity interval versus 
relative velocity. Since Ovi A 1037 is heavily blended we use only 
the O vi/11031 line. For Civ, we have used the oscillator strength 
weighted mean of the column densities per unit velocity measured 
from both absorption lines of the doublet. For clarity, we have mul- 
tiplied the Civ apparent column density profile by a factor of 10. 
Vertical dashed and dotted lines show the positions of peaks in the 
optical depth of C iv and O vi respectively. It is apparent that the 
O vi peaks are shifted compared to the C iv ones. 

In the lower panel of Fig. [5] we plot the ratio of O vi to C iv 
apparent column densities per unit velocity against the relative ve- 
locity and find that the ratio varies between 10 and 20 through the 
Civ absorption profile. The fact that the Ovi absorption profile is 
broader suggests the existence of gas outside the C iv profile with 
O vi to C iv column density ratio higher than 20. 

The component at — \-4.2 km s~'has virtually no detectable H 1 
absorption associated. From the Lyman-ff line we derive an upper 
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Figure 4. Absorption profiles are shown on a velocity scale for regions with low r\ values (see Fig. 1): D, F (upper-left and right hand-side panels) and H, 
N (lower left and right hand-side panels). C iv and O vi are clearly detected in regions D and F. In these cases we also overplot the best fitted Voigt profile 
model and indicate individual components with vertical tick marks. For regions H and N, there are only tentative O vi coincidences. Vertical dashed lines in the 
bottom panels mark the locations of tentative O vi doublet components. The vertical dashed lines in the upper-left panel delineate the region of O vi absorption. 
Regions marked by 'x' are Lyman-a contamination. 



limit of log N(H i) = 12.80 suggesting that metallicity is probably 
high in this component. Indeed, given the low b value of the com- 
ponent, it is probable that the gas is photo-ionized in which case 
metallicity has to be close to solar. For the other three components 
that coincide with a strong H i absorption and it is impossible to 
quantify the amount of H i absorption associated with them indi- 
vidually such that useful metallicity limits can be established. 

We have seen before (Fig. O that the x 1 curve correspond- 
ing to the fit of the He n absorption shows a marked minimum for 
17= 12 in the case of turbulent broadening (i.e £ = fcH elI /^Hi = 1) 



and no minimum for pure thermal broadening = 0.5). In Fig. [6] 
we show in the left hand-side row the simulated He 11 profiles for 
?7~ 130 and £ = 0.5 (solid curve) and 77 ~ 12 and £ = 1.0 (dashed 
curve). Remember that for these fits we have used the minimum 
number of Voigt profile components without any constraint from 
the O vi profile. It is apparent that the red solid Hen profile (ob- 
tained assuming pure thermal broadening) is missing several pixels 
in the red wing of the region of interest around km/s. This is be- 
cause the b value of the corresponding component (fixed by the 
H 1 profile) must be much larger to reach this position. If we now 
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Table 2. Component Detail of the System at z a bs = 2.6346 
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v (kms ') 



Ion 



H i lines used 



> (km s 



log (N in cm 2 ) 
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Figure 5. Top panel : Apparent column density profiles of O vi (in blue) and 
C iv (in red) per unit velocity interval versus relative velocity for the system 
z a bs = 2.6498. For display purpose, C iv apparent column densities have 
been multiplied by a factor of 10. The mean errors in each pixel for Ovi 
(in blue errorbar) and C iv (in red errorbar) are shown in the left. Bottom 
panel : Apparent column density ratio through the profile. The errorbar in 
the left shows corresponding mean error in each pixel. Vertical dotted lines 
show the centroids of Voigt components. The shift between the O vi and 
C iv centroids is apparent. 



Table 1. Results of multiple component Voigt profile fitting to the z a b s = 
2.6498 system. The parameters of the H i are obtained by keeping the com- 
ponent structure as seen in Ovi. 



v (km s 1 ) 


Ion 


H i lines used 


b(km s _1 ) 


log (N in cm 1 ) 


-214.6 ± 1.0 


Hi 


Ly-or 




16.3 ± 1.5 


12.59 ± 0.03 


-150.8 ± 0.3 


Hi 


Ly-or, 


Ly-/? 


23.6 ±0.3 


14.63 ± 0.01 


-76.2 ± 0.0 


Hi 


Ly-or, 


Ly-/? 


37.6 ± 0.9 


14.91 ± 0.01 


-48.6 ± 0.0 


Hi 


Ly-a, 


Ly-/? 


23.1 ±0.7 


14.68 ± 0.02 


-20.2 ± 0.0 


Hi 


hy-a, 


Ly-/? 


26.4 ± 1.1 


13.53 ± 0.06 


-76.2 ± 3.3 


Ovi 






38.6 ±3.9 


13.56 ± 0.05 


-58.8 ± 1.0 


Civ 






14.5 ± 1.5 


12.30 ± 0.03 


-48.6 ± 0.3 


Ovi 






9.4 ± 0.9 


13.16 ±0.06 


-24.1 ±0.6 


Civ 






5.9 ± 1.0 


12.06 ± 0.04 


-20.2 ± 0.2 


Ovi 






10.9 ± 0.4 


13.55 ± 0.02 


+ 1.9 ±0.4 


Civ 






4.8 ±0.7 


12.16 ±0.03 


+4.2 ±0.1 


Ovi 






6.5 ±0.2 


13.49 ± 0.01 



add the constraint that H i should be associated with the three O vi 
components, we can reproduce this profile better. Indeed, because 
of the extra component at v ~ -20 km/s, pure thermal broadened 
Hen profile with higher 77 (~ 100) gives an equally good fit (right 
panel of Fig. 6). 

It seems therefore that if we add a He n component at the po- 
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sition of the redder O vi component, any value of rj between ~15 
and 100 is acceptable. Thus it seems that the possible presence of 
unresolved narrow H 1 components could be one of causes of low j] 
measurements. It is a fact however that the main H 1 components 
have large b values, corresponding to temperatures in excess of 
10 5 K. Therefore it is not impossible that part of the gas is at high 
temperature. 



4.2 System at z abs = 2.6346 

Absorption profiles from this system are shown on a velocity scale 
in Fig. [4] Unlike in the previous system the velocity range of metal 
lines falls well within the Lyman-a profile. The Ovi/11037 line is 
blended with Lyman-y at z = 2.8781 and Lyman-/? at z = 2.6765. 
Because of this contamination we use the well measured redshifts 
of Civ components to fit the Ovi doublet. The contributions of 
the contaminating lines are self-consistently computed using other 
available transitions. In addition to the C iv counterparts, we need 
two components in the red part of the profile to fit the O vi dou- 
blet where there is no C iv absorption. H 1 Lyman-o-, Lyman-/?, and 
Lyman-y lines have been fitted simultaneously imposing compo- 
nents at the redshifts of five O vi components. Two extra compo- 
nents are required in the blue ( 100 km/s) to cover the total Hi 

absorption. The details of the fit results are given in Table[2] 

As in the previous system, for the components with both C iv 
and O vi the O vi ^-parameters are larger than the C iv ones and the 
column density ratio of Ovi to Civ is as high as ~ 15. The Ovi 
^-parameters correspond to upper limits on the kinetic temperature 
of 4 x 10 5 , 3 x 10 s and 9 x 10 4 K respectively for the components 
at -16.4, -0.6 and +36.2 km s -1 . In the components where we find 
only O vi the ratio of O vi to C iv column densities can be higher 
than 20. These components have broad O vi lines with ^-parameters 
corresponding to upper limits of 6 x 10 5 and 10 6 K respectively for 
the components at +48 and +101 km s~ x . The corresponding H 1 
components also have high b values allowing for high temperature 
(~ 10 5 K) in the gas associated with these two components. All 
this suggests a multiphase structure in this absorbing gas with the 
possible existence of a hot phase contributing to most of the O vi 
absorption. Indeed, the O vi profile is suggestively broad. 
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Figure 6. Fits of H i and He n absorption in the z-abs, = 2.6498 system. In all panels dashed (blue) and solid (red) curves are for turbulent, f = fc(HeII)/fo(HI) = 1, 
and thermal, ^ = 0.5, broadening cases respectively. Left hand-side column: Fit using the minimum number (four) of components required to fit the H i 
absorption. The dashed and solid curves in the top panel are the simulated Hen profiles with n ~ \2 and 130 respectively. Only pixels between the two vertical 
dot-dashed lines are used to derive f; by x 2 minimization. Middle column : Fit with five components imposing H i components at the position of the O vi 
components. The dashed and solid curves in the top panel are the simulated Hen profiles with i; ~ 12 and 100 respectively. Note that two components around 
+90 km/s are fitted together to take into account their contribution to the He n absorption. For these two components we take ;/ = 65 and pure turbulent 
broadening derived from the analysis of region F' (see Fig. [2). Right hand-side column: x 2 plot for fits shown in the middle row. Vertical ticks in middle and 
left panels indicate the positions of individual Voigt components. 



We fitted the H i and He n profiles in the two extreme cases 
of pure turbulent and pure thermal broadening, considering both 
components from the fit of the H i profile only and from the fit of 
metal lines. Results are given in Fig|7]and Table 2. We notice from 
right panel of Fig [7] that even when we tie the H i components to 
O vi components the best fitted x 2 is obtained for the pure turbu- 
lent case with low 77. However, reality probably corresponds to an 
intermediate case with <f between 0.5 and 1. In the bottom panel of 
Fig. [8] we plot the minimum^ 2 value obtained for different values 
of £ . Even though the best fitted x 2 value is obtained for £ = 1 the 
curve is flat and the lcr range is if > 0.6. As can be seen in the top 
panel of the figure, this can accommodate a wide range of 77. There- 
fore, in this system also high 77 values are acceptable although H 1 
and O vi absorption profiles are broad and highly suggestive of a 
gas with temperature higher than the typical photo-ionization tem- 
perature (i.e few 10 4 K). 



5 MODELS 

Given the particularities of the systems singled out by the pres- 
ence of O vi absorption, possible low 77 values and high O vi/C iv 
ratios, we have constructed models to test the different mecha- 
nisms that could induce such properties. It is well known that 
photo-ionization by a power-law spectrum with appropriate slope 
can yield low 77 values. This would require the presence of local 
sources of hard phot ons (see IShullet al.ll2004h . Observations by 
IWorseck et al. ( 2007) seem to show however that there is no QSO 
present in the vicinity of the two absorbers considered in the previ- 
ous Section. While this observation does not rule out a QSO emis- 
sion highly beamed perpendicular to the line of sight or a short lived 
QSO emission in the vicinity of the absorbers, we explore alternate 
explanations for low 77 in the O vi absorbers. Therefore, in the fol- 



lowing we present the results of models of a hot gas embedded in 
the meta-galactic UV background. 

We use the photo-ionization code Cloudy (v07.02; 
iFerland et alj 0998)) to derive the ionization structure in a 
gas with fixed temperature (therefore not controlled by photo- 
ionization). This will make it possible to discuss at the same 
time both extreme situations (collisional ionization and photo- 
ionization) but also the intermediate situation of high-temperature 
gas with a contribution of photo-ionization. For comparison, we 
also show results from the model in which the temperature is the 
consequence of thermal equilibrium under photo-ionization. The 
calculations are made in the optically thin case. We use the Haardt 
and Madau (2005) background spectrum dominated by QSOs. 
We assume relative solar abundances and [C/H] = -1.0. In the 
top panel of Fig. [9] we plot the variation of the O vi to C iv ratio 
with hydrogen density. The solid black line is the result of model 
calculations where temperature is calculated by CLOUDY assum- 
ing photo-ionization equilibrium. Other lines are for temperatures 
in the range 5xl0 4 -5xl0 5 K. It is to be remembered that when 
pure collisional excitation is considered the fraction of He n is 
maximum whe n 4.5 < log T(K) < 4.9 and in the case of O vi it is 
T ~ 3x10 s K dGnat & Sternberg||2007l) . At low temperature (say 
T < 5 x 10 4 K) the ionization is dominated by photo-ionization. As 
expected the transition between photo-ionization dominated and 
collisional ionization dominated regimes happens at T ~ 10 5 K. 

The horizontal dotted lines show the range of observed O vi 
to Civ column density ratios (between 10 and 20) seen in the 
Civ components of the two systems discussed above. This range 
is well reproduced by models with T< 10 5 K for a typical den- 
sity of 10~ 4 cm" 3 . However the higher O vi to C iv ratio inferred in 
the velocity range (or Voigt profile components) where only O vi 
is detected needs either low density (and low temperature) photo- 
ionized gas or high density (i.e > 10~ 3 ctrr 3 ) hot gas (T> 10 5 K) 
where collisions begin to play a role. Interestingly such high tem- 
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Figure 7. Fits of H i and He n absorption in the z-jbs = 2.6346 system. In all panels dashed (blue) and solid (red) curves are for turbulent, f = fc(HeII)/fo(HI) = 1, 
and thermal, f = 0.5, broadening cases respectively. Left hand-side column: Fit using the minimum number (five) of components required to fit the H i 
absorption. The dashed and solid curves in the top panel are the simulated Hen profiles with n ~ \2 and 160 respectively. Only pixels between the two vertical 
dot-dashed lines are used to derive n by x 2 minimization. Middle column: Fit with seven components imposing H i components at the position of the O vi 
components. The dashed and solid curves in the top panel are the simulated Hen profiles with n ~ 7 and 140 respectively. Right hand-side column: x 2 plot for 
fits shown in the middle row. Vertical ticks in middle and left panels indicate the positions of individual Voigt components. 



Region D | z abs = 2.6346 




Figure 8. System at z a bs= 2.6346. Bottom: Minimum^ 2 as a function of f 
in the case of H i fitted with minimum number of components. Top: Best 
fitted value of r\ for different values of f . 



peratures are not ruled out by the fo-parameters of O vi components 
(see discussions in the previous Section). 

In the bottom panel of Fig. [9] we plot r] as predicted by the 
models versus the hydrogen density. It is apparent that low rj values 
(i.e < 60) are only possible for T > 10 5 K. Available data on H i 
and He n profiles allow for the existence of such hot gas that would 
also produce the component with high O vi to C iv column density 



ratio (i.e N(0 vi)/N{C iv)>20). It is apparent from Fig. 4 that the 
absorption profiles indicate higher Doppler parameters going from 
C iv to O vi to H i. This has already been noted for C iv and O vi 
by Fox et al. (2007) and interpreted as the existence of a hot phase. 
We note that the b values measured for the strongest H i compo- 
nents in the two systems (38.6 km/s at z abs = 2.6498 and 46.6 and 
69.5 km s -1 at z abs = 2.6346, see Tables 1 and 2) are consistent with 
a temperature, T > 10 5 K and it is apparent from the absorption 
profiles that larger b values could be accommodated. 

If the low t] values were to be confirmed, we would favor a 
scenario where the absorbing gas is a multiphase medium in which 
photo-ionized gas components coexist with a wide range of den- 
sity and temperature. While most of the metal absorption traced by 
C rv comes from relatively cold (i.e T< 10 5 K) gas part of O vi and 
predominant contributions of H i and He n could be due to a hot 
phase (T > 10 5 K). There is evidence for the existence of mul- 
tiphase media in the low-z O vi absorbers l lTripp et aT]|200£ ) and 
O vi absorption associated with high-z DLAs jFox et alj|2007 ). 



6 CONCLUSIONS 

We have reanalyzed the line of sight towards the high redshift 
fom = 2.885) bright quasar QSO HE 2347-4342 and measured the 
parameter rj = N(He n)/N(H i) in the Lyman-a forest using Voigt- 
profile fitting of the H i transitions in the Lyman series. As in pre- 
vious studies, we find that ^ > 50 in most of the Lyman-a forest 
except in four regions where it is much smaller (77 ~ 10 — 20). 

We detect O vi absorption associated with two of these re- 
gions (at z a bs = 2.6346 and 2.6498). The corresponding wavelength 
ranges for the two other regions are too blended to reach any firm 
conclusion on the presence of associated O vi absorption. We ob- 
serve that the z abs = 2.6346 system is a usual system with the metals 
located at the center of the H 1 profile whereas the z a bs = 2.6498 sys- 
tem has the metals displaced in the red wing of the H 1 absorption 
but moreover, with the C iv profile systematically shifted compared 
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Figure 9. Top panel: The O vi to C iv column density ratio versus hydro- 
gen density. The different curves correspond to the results of optically 
thin CLOUDY models with constant temperature except the solid black 
curve which is for thermal equilibrium. The gas is assumed of metallicity 
[C/H] = —1.0 & solar relative abundances and exposed to QSO dominated 
HM05 ionizing flux. The horizontal dotted and dashed lines show the range 
of observed column density ratios for the system at Zubs = 2.6498 and 2.6346 
respectively. The vertical dotted and dashed lines indicate the corresponding 
allowed range in hydrogen density for photo-ionization equilibrium. 

to O vi. Doppler parameters of the well-defined C iv components 
rule out the fact that the associated gas is hot and favor the idea that 
this gas is photo-ionized. We show that if we constrain the fit of 
the H i and/or He n absorption profiles with the presence of metal 
components, we can accommodate r] values in the range 15-100 in 
these systems assuming broadening is intermediate between pure 
thermal and pure collisional. 

We construct constant density photo-ionized models and show 
that while simple photo-ionization models reproduce the observed 
./V(0 vi)/iV(C iv) ratio for a range of density, they fail to produce 



low rj values. On the contrary, models with high temperature (i.e 
T > 10 5 K) can produce low values of rj. The Doppler parameters 
of the strongest H i components are consistent with such a tem- 
perature. In addition, the higher b values observed for Ovi com- 
pared to C iv and the existence of O vi alone components suggest a 
multiphase nature of the absorbing region. Therefore, if low rj val- 
ues were to be confirmed, we would favor a multi-phase model in 
which most of the gas in the regions of low rj is at high temperature 
(>10 5 K) but the metals and in particular C iv are located in lower 
temperature photo-ionized and probably transient regions. As the 
high temperature gas can not be produced by photo-ionization, we 
expect the O vi systems with low rj to be associated with galaxies. 
Therefore, deep search for Lyman break galaxies at these redshifts 
may be interesting to perform in these fields. 
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